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ABSTRACT 

Studying 3 samples of early- type galaxies, which include approximately 8800 
galaxies and cover a relatively ample magnitude range (< AM > ~ 5.5 mag), we 
find that the coefficients as well as the intrinsic dispersion of the Fundamental Plane 
depend on the width and brightness of the magnitude range within which the galax- 
ies are distributed. We analyse this dependence and the results show that it is due 
to the fact that the distribution of galaxies in the space defined by the variables 
log(re), < /I >e, log((To) dcpcuds on the luminosity. 

Key words: Galaxies: elliptical and lenticular, fundamental parameters. 



1 INTRODUCTION 

The Fundamental Plane (FP) for ordinary early-type 
galaxies (ETGs) is a correlation between the following 
observable parameters: logarithm of the effective radius 
(log (re)), effective mean surface brightness (< /i >e) 
and logarithm of the cen t ral velocity dispersi on (log ctq) 
l|Diorgovsky fc Davis 19871 : iDressler et al. 19871 ). This cor- 
relation is usually expressed as follows: 



a + (3 log(re). 



(2) 



Different studies have demonstrated that ETGs 
inside clusters define the KR with an intrinsic dis- 
pers ion of approximate ly 0.4 ma g in < fi >e 



(.ffamabe Sz Kormendv 1987i : 
ISandage k. Peremulter 199ll : 



Ifoessel et al. 19871 : 
ISandage fc Lubin 200lf ). 



log (re) = a log ((To) + b < a* >e + c. 



(1) 



where a, b and c represent scale factors. These scale fac- 
tors contain important information about ETGs and might 
depend, among other things, on environment and wave- 
length. 

The FP relation (equation 1) is a consequence of the 
dynamical equilibrium condition (virial theorem) and of 
the regular behaviour of the mass-luminosity ratio as well 
as the structure along the whole range of luminosities for 
the ETGs. Due to the small intrinsic dispersion (~ 0.1 dex 
in re and (tq, and ~ 0.1 in < /i >e), the FP is considered a 
powerful tool for studying galactic formation and evolution 
as well as useful for making est i mations of large scale 
distances (iKiaergaard et al. 1993: Jorgensen et al. 19961 : 



IJorgensen" et al. 19991 : iKelson et al. 19971 ) 

An important projection of the FP is the correlation 
between log(re ) and < fj, >e, al so known as the Kormendy 
relation (KR) (jKormendv 19771 '): 
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Recent studies have demonstrated that ETGs contained in 
less dense environments and isolated ETGs also follow the 
KR with t he same coeffi cients and dispers ion as do ET Gs 
in clusters l|Reda et al. 2004 : Nigoche-Netro et al. 2007h . 

On the other hand, Nigoche-Netro (2007) and Nigoche- 
Netro et al. (2007, 2008) found, using observational data 
and numerical simulations, that both the coefficients and 
the intrinsic dispersion for the KR depend on the width 
of the magnitude range and the brightness of ETGs 
within the magnitude range. This dependence is caused 
by a "geometricai effect" due to the fact that the dis- 
tribution of ETGs on the log(re) - < fj, >e plane de- 
pends on luminositv ijVarela 20041 : iD'Onofrio et al. 20061 : 
iNigoche-Netro et al. 2007i ) and that the geometric shape of 
this distribution is not symmetrical. In other words, the ge- 
ometric shape of the distribution of ETGs on the log(re) - 
< >e plane changes systematically as brighter ETGs are 
progressively considered, and so the values of the KR coef- 
ficients change too, because the fitting of a straight line to 
a set of data does not give the same result for slope and in- 
tercept values for data distributed with a rectangular shape 
as for data distributed with a triangular shape or, for that 
matter, with another geometrical shape. 

The previously mentioned results, together with the fact 
that the KR is a projection of the FP, suggested that a study 
of the behaviour of the coefficients and intrinsic dispersion 
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of the FP with respect to the magnitude range would prove 
interesting. Several papers in the literature (Bernardi et al. 
2003c; Jorgensen et al. 1996) undertake, somewhat differ- 
ently, this type of investigation. These works have inves- 
tigated, among other things, the residues of the FP with 
respect to magnitude. The results they find indicate that 
there is no correlation between the residues of the FP and 
magnitude. However, we consider that it is necessary to per- 
form an investigation similar to that presented for the KR 
in Nigoche-Netro et al. (2008) for the case of the FP, given 
that there is no published work in the astronomical litera- 
ture where this approach is taken. 

This study proves to be relevant for investigations that 
perform comparisons of different galaxy samples such as 
comparisons of samples for which dependence on the en- 
vironment, on the redshift or on the wavelength are sought 
for. If differences are found and the dependence of the coeffi- 
cients of the FP relation on the magnitude range is not taken 
into consideration, the differences found may be completely 
misinterpreted. 

In this paper we compile 3 samples of ETGs which in- 
clude a total of approximately 8800 galaxies and cover a 
relatively ample magnitude range (< AM > ~ 5.5 mag). It 
is important to mention that in this work we perform the 
analysis of the coefficient values and of the intrinsic disper- 
sion of the FP with respect to: 1) the width of the magnitude 
range (increasing magnitude intervals), that is, we consider 
the faintest galaxies in each sample and we progressively in- 
crease the width of the magnitude interval by inclusion of 
the brighter galaxies, and 2) the brightness of the magni- 
tude range (narrow magnitude intervals), that is, we con- 
sider galaxy samples in progressively brighter fixed-width 
magnitude intervals. On the other hand, in order to avoid 
biases in the coefficients, we use the linear regression method 
known as Measurement errors an d Intrinsic Scatter Thre e 
dimensional bisector (MISTsis) l|La Barbera et al. 2000h . 
This method takes into consideration the errors associated 
to each one of the variables and minimises the effects of the 
intrinsic dispersions of each galaxy sample under study. 

In Section 2 we introduce the different galaxy samples. 
In Section 3 we present the fitting method used for calculat- 
ing the parameters of the FP and the results from these fits. 
In this section we also discuss an analysis of the behaviour 
of the coefficients a, b and c and of the intrinsic dispersion 
for the FP with respect to the absolute magnitude range. 
Finally, in Section 4, we present our conclusions. 



2 THE SAMPLES 

We use a sample of 8666 ETG s from the Sloan Digit al Sky 
ourvey (total SDSS sample) l|Bernardi et al. 2003al '). This 
sample has surface photometry information in filters g*, r*, 
i* and z* (absolute magnitude range — 19 ^ Mr* > —25 
and its equivalent in other filters). We also use a sample 
of 116 ETGs for the Coma Cluster l|Milvang- Jensen 19971 : 
IJorgensen et al. 19991 ) in the filter Gunn r (-20 ^ Mcr > 
-24.9) and a 44 ETGs sample in th e filter Gunn r (-19 ^ 
Mar > -24) from the Hydra cluster (|Milvang- Jensen 19971 ). 

All the samples are redshift-homogeneous (the galax- 
ies are contained within a narrow redshift interval), ex- 
cept for the total SDSS sample which cover a relatively 



ample redshift interv al (0.01 z ^ 0.3 ). This sample 
is magnitude- limited (jBernardi et al. 2003bl ). that is to say 
that in small volumes and for small values of z, it does 
not contain very bright galaxies, whereas in large volumes 
and for large z values, it does not contain faint galaxies. 
Besides, within large volumes there could be evolution ef- 
fects of the parameters of the galaxies. So, in order to have 
a representative sample of the universe in a given volume 
without any evolution effects it is important to consider 
narrow redshift intervals. Bernardi et al. 2003b recommend 
Az = 0.04. This value comes from the sizes of the largest 
structures in the universe seen in nu merical simulations o f 
the cold dark matter family of models (jColberg et al. 20001 ). 
It is also well known that the SDSS photometry underesti- 
mates the luminosity of the brightest objects in crowded 
fields (|Bernardi et al. 2007al : iBernardT et al. 2007bl ). In ad- 
dition, the spectroscopic reduction of the SDSS data under- 
estimates the valu es of the small central v elocity dispersions 
((TO s; 150 km/s) iBernardi et al. 2007bh . 

To probe the possible evolution effects and avoid the 
photometric and spectroscopic biases we have built a sub- 
sample from the total SDSS in the redshift interval 0.04 gC 
z ^ 0.08 and central velocity dispersion ctq > 150 km/s. 
Given that the determination of the distances to the galax- 
ies relies entirely on z, we chose the lower limit to keep the 
inffuence of the peculiar velocities on the measured z below 
10%. This subsample has 933 galaxies in each filter and cover 
a magnitude range < AM > ~ 4 mag (-19 > Mg* > -22.8 
and its equivalent in other filters). From now on, this sub- 
sample will be referred to as homogeneous sample from the 
SDSS. 

In Table 1 we present relevant information (number 
of galaxies, magnitude range and redshift) for the samples 
of galaxies we use in this paper. The magnitude range in- 
formation is given in relation to the different filters used 
(from the literature) and also the approximate range for 
the B- magnitude (calculated by us). The transformation 
to the B filter was accomplish ed by use of the fol lowing 
equation: B - Gunn r — 1.15 (|Milvang- Jensen 19971 ). The 
filters g*, r*, i* and z* correspond approximately to the 
Johnson- Morgan-Cou sins filters B, V, Rc and Ic respectively 
(|Fukugita et al. 19961 ). 

All the samples consider the photometric and spectro- 
scopic parameters log(re), < ^ >e and log(CTo) corrected for 
different biases (galactic extinction, K correction, cosmolog- 
ical dimming and aperture correction). These parameters 
as well as their uncertainties were taken directly from the 
different papers cited above. It is worth stressing that the 
photometric parameters were obtained using de Vaucouleurs 
r^^^ fits. These fits represent a good approximation to the 
Sersic profile- fits because all galaxies used here are bright 
{Mb s; -18) (|Prugniel fc Siemen 199?! ). To stress this point 
we should mention that Nigoche-Netro et al. (2008) show 
that the de Vaucouleurs approximation does not produce 
appreciable biases in the behaviour of the KR parameters 
with respect to the magnitude range. 

Finally, an important characteristic of SO galaxies is 
that, in general, the structural properties of their bulges 
show approximately the same homogeneity as those of E 
galaxies. Due to this fact, and because the FP uses effective 
parameters, we use the photometric parameters from the 
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Table 1. Name, number of galaxies, magnitude range (according to the papers where the samples were taken from), a 
pproximate magnitude range in the B-Filter (calculated by us) and rcdshift of the different galaxy samples compiled in this 
work. 



Sample 




N 


Magnitude 


rant 




Approximate magnitude range in the B-filter 


z 


Total SDSS (g* filter) 




8666 


-18.0 ^ A/g* 


> - 


24.1 


-17.5 ^Mb> 


-23.6 


< 0.3 


Total SDSS (r* filter) 




8666 


-18.6 > A[r* 


> - 


-24.7 


-17.5 ^Mb> 


-23.6 


< 0.3 


Total SDSS (i* filter) 




8666 


-19.0 ^ Af„ 


> - 


25.1 


-17.5 ^Mb> 


-23.6 


< 0.3 


Total SDSS (z* filter) 




8666 


-19.3 > M^* 


> - 


25.3 


-17.5 ^Mb> 


-23.6 


^ 0.3 


Coma cluster (Gunn r 


filter) 


116 


-20.4 ^ Mcr > - 


^24.9 


-19.3 ^ Mb > 


-23.8 


0.024 


Hydra cluster (Gunn r 


filter) 


44 


-19.0 ^ Mcr 


> - 


-24.6 


-17.9 ^ Mb > 


-23.5 


0.014 



whole galaxy in the case of bright Es, and from the bulge in 
the case of bright SOs. 



3 THE FP RELATION FOR THE DIFFERENT 
SAMPLES OF ETGS 

3.1 Calculation of the FP relation 

The estimation of the FP coefficients may be severely af- 
fected by the fitting method used and by the independent 
variable adopted for the fit. The biases may be larger if 
there are measurement errors in the variable, if these errors 
are correlated or if there is intrinsic dispersion. The Measure- 
ment errors and Intr insic Scatter Three dime nsional bisector 
method (MISTsis) (|La Barbera et al. 200o l) is a statistical 
model that takes into consideration the different sources of 
bias mentioned above. Here we use this method to calculate 
the FP coefficients. This fitting method is a 3 dimensional 
generalisation of the Bivariate Correlated Errors and Intrin- 
sic Scatter bisector [BCESbis) method described by Isobe 
et al. (1990) and Akritas & Bershady (1996). The coeffi- 
cients for the edge-on FP are calculated using the BCES Bis 
method. 

From the photometric and spectroscopic parameters for 
the different galaxy samples we calculate the FP coefficients 
in increasing and narrow magnitude intervals (see Section 
1), the errors of these coefficients and the total intrinsic dis- 
persion (cr) of the FP relation (subtracting in quadrature 
from the intrinsic dispersion in < /x >e the residues disper- 
sion due to the measurement errors of < fj, >e, log(re) and 
log(cro)). The definition of the total intrinsic dispersion of 
the FP is a generalisation in three dimensions of the defi- 
nition of the total intrinsic dispersion of the KR given by 
La Barbera et al. (2003). According to La Barbera et al. 
(2003), for the calculation of the total intrinsic dispersion it 
is necessary to have the measurement errors of each one of 
the variables involved. These errors were taken from the lit- 
erature (see Section 2.1), while the theoretical errors in the 
FP coefficients were calculated by us following La Barbera 
et al. (2000) in la intervals. 

In Tables 2 and 4 we present the results for the values of 
the coefficients of the FP for the different samples of galaxies 
and for different magnitude intervals. 



3.2 Behaviour of the FP coefficients with respect 
to the width of the magnitude range 

As may be seen from the results given in Table 2 for the 
different galaxy samples, the intrinsic dispersion increases 
when a growing number of brighter galaxies is included in 
the samples (upper magnitude cut-off), that is to say when 
we use increasing magnitude intervals. We also observe that 
the a, b and c coefficients increase systematically as brighter 
galaxies are included (see Fig. 1) and that these changes 
are larger than the associated errors for most of the cases 
(differences between a coefficients may be as large as 33%). 
On the other hand, if we consider the brightest galaxies in 
each sample and we progressively increase the width of the 
magnitude interval by inclusion of the fainter galaxies (lower 
magnitude cut-off), the behaviour of the intrinsic dispersion 
and FP coefficients is similar to the one described above, 
that is, the FP parameters increase systematically when we 
increase the width of the magnitude interval. 

To investigate whether the changes in the coeffi- 
cients are significant or whether they are only the prod- 
uct of statistical ffuctuations, it is necessary to apply 
non-parametric tests to the data for the different galaxy 
samples. Among the non-parametr ic tests appropriate for 
our study we find the r un test (|Bendat k, Piersol 19661 : 
iNigoche-Netro et al. 200l ). 

In Table 3 we show the results of the application of 
the run test method to the data for coefficients a, b and 
c. We take as a null hypothesis that there is no underlying 
trend in the data. The percentages given in Table 3 refer 
to the confidence level with which we may reject the null 
hypothesis. From Table 3 we may see that, on average, the 
null hypothesis may be rejected with a confidence level of 
96%. This implies that there are strong reasons to affirm 
that there is an underlying trend in the values of each one 
of the coefficients of the FP. 

In appendix A we present an analysis of the robust- 
ness of the run test using the bootstrap methodology 
(|Efron 19791 ). This analysis includes measurement errors as- 
sociated with the data. The conclusions of this analysis indi- 
cate that the run test is very efficient in finding underlying 
trends in sets of data even if the errors associated are quite 
large. 

Here, it is important to mention that when we compare 
the behaviour of the FP coefficients for the total SDSS sam- 
ple and the homogeneous SDSS sample (Fig. 1) we note that 
this behaviour is similar for both samples, meaning that the 
coefficients for the homogeneous sample increase systemat- 
ically as we increase the width of the magnitude range. It 
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Table 2. FP coeiRcients for the different samples of galaxies in increasing intervals of magnitude 
(upper magnitude cut-off). MI is the absolute magnitude interval within which the galaxies are 
distributed, AM is the width of the magnitude interval, N is the number of galaxies in the 
magnitude interval, a, 6 and c axe the FP coefficients (obtained using the MIST^jj method) and 
a is the total intrinsic dispersion of the FP. 
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Figure 1. Behaviour of the FP coefficients with respect to the width of the magnitude range (upper magnitude cut-off). Circles represent 
the total SDSS sample (g* filter). Diamonds represent the homogeneous SDSS sample (g* filter). 



Table 3. Run test for the evaluation of the 
FP coefficients (increasing magnitude intervals) 
from the different samples of galaxies. The null 
hypothesis establishes that there is no underly- 
ing trend in the data, the percentages refer to 
the confidence level with which we can reject the 
null hypothesis. 



Sample 


a 


b 


c 


Total SDSS (g* filter) 


98% 


98% 


98% 


Total SDSS (r* filter) 


98% 


98% 


98% 


Total SDSS (i* filter) 


98% 


98% 


98% 


Total SDSS (z* filter) 


98% 


98% 


98% 


Coma (Gunn r filter) 


93% 


93% 


93% 


Hydra (Gunn r filter) 


93% 


93% 


93% 



is also clear that the values of these coefRcients result in 
larger values than those for the total sample. However, it 
is not possible to say that this diference may be caused by 
evolution, because the biases in the structural parameters, 
and the cut in the velocity dispersion might also have an 
effect on the values of the coefficients (see Section 3.5 for 
full details). 



3.3 Behaviour of the FP coefRcients with respect 
to the brightness of the magnitude range 

In order to analyse the behaviour of the FP coefficients with 
respect to the brightness of the magnitude range and study 
the distribution of the galaxies in the space defined by the 
log(re), < ^ >e, log(cr) variables, we may utilise the FP edge- 
on graphs. Figure 2 shows the results for the total SDSS sam- 
ple in the g* filter, where each symbol and colour represent 
a one-magnitude wide interval and the thick continuous line 
represents the fit to all the data. Table 4 shows the values 
of the FP coefficients. These values were obtained from the 
fits {BCESsis) made to the different magnitude intervals, 
however in order to analyse the effects that the different in- 
dependent variables may have on the results of the fits, we 
present the results using as independent variables log(re) 
and < fi >e- 

In Figure 2 and Table 4 we can see that different mag- 
nitude galaxies are progressively displaced from the lower 
part of the diagram (faint galaxies) towards the upper part 
(bright galaxies) and that the values of the coefficients and 
of the dispersion of the FP change systematically. These 



changes are larger than the associated errors for most of 
the cases (application of the run test confirms that there 
exists an underlying trend in the data for each one of the 
coefficients) . It is also clear that the fit to the complete sam- 
ple shows significant differences with the fits performed to 
galaxies of different luminosity. For example, the differences 
between coefficients a, b and c for the middle part of the 
galaxy distribution (-21.0 M >-22.0) and the a, b and c 
coefficients for the total sample (-18.5 ^ M >-24.0) is ap- 
proximately 20%. From what we have just said it can be 
concluded that different luminosity galaxies define consec- 
utive planes with different values for the coefficients and 
the dispersion. This implies that the distribution of galax- 
ies in the space defined by variables log(re), < ^ >e,\og{a) 
changes systematically as brighter galaxies are considered 
(see Section 3.5 for full details). 

What we have just described goes in the same direc- 
tion as what Bender et al. (1992) reported which indicates 
that the distribution of different stellar systems (elliptical, 
spiral bulges, compact elliptical and dwarf ellipticals) in the 
parameter space log(re), < >e, \og{(j) depends on luminos- 
ity. According to Bender et al. (1992), this means that the 
different systems define planes on this space with possible 
small, parallel offsets between systems. 

An alternative way for studying the behaviour of the FP 
as a function of the magnitude range, is to analyse the values 
and distribution of the residues, that is to say, to analyse the 
difference between the predicted values from the best fit to 
the FP data and the real values. However, given that the 
values of the coefficients of the FP for the entire sample are 
different from the values of the coefficients of the samples in 
different magnitude ranges, we shall analyse the behaviour 
of the residues taking the coefficients both for the entire 
sample (-18.5 ^ M >-24.0) as well as for the middle part 
of the distribution (-21.0 ^ M >-22.0 ). In Figures 3 and 4 
we show the orthogonal residues of the FP (g* filter) with 
respect to the log(re), < >e, log((j), M and Xpp variables. 
Where Xpp represents the distance along the major axis of 
the FP (see Bernardi et al. 2003c). 

As can be seen from Figure 3, when the residues are 
obtained considering the coefficients of the FP for the com- 
plete sample, the distribution of the residues with respect 
to the different variables which we have used show no corre- 
lation except, perhaps, for the case of the log((T) and Xpp 
variable. However, when we consider the FP coefficients for 
the middle part of the galaxy distribution (-21.0 ^ M >- 
22.0) (Figure 4) we find that the residue distribution with 
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log(a ) + 0.2 <^> 



Figure 2. Distribution of the galaxies in the complete SDSS sample (g* filter) on the edge-on plane. Each symbol and colour represent 
a one-magnitude wide interval. The continuous line is the fit (BCESsia) obtained for all the data in the sample. 



Table 4. FP coefficients for the total SDSS sample (g* filter) in narrow magnitude intervals. MI 
is the absolute magnitude interval within as the galaxies are distributed, AAI is the width of the 
magnitude interval, N is the number of galaxies in the magnitude interval, a, b and c are the FP 
coefficients and a is the intrinsic dispersion of the FP. 



MI 


AM 


N 


a 


b 


c 




BCESsis 


-18.5 ^ M >-19.0 


0.5 


26 


1.145±0.050 


0.229±0.002 


-6.874±0.315 


0.240 


-19.0 ^ M >-20.0 


1.0 


492 


1.157±0.029 


0.231±0.001 


-6.862±0.185 


0.294 


-20.0 ^ M >-21.0 


1.0 


2341 


1.126±0.012 


0.225±0.001 


-6.603±0.075 


0.289 


-21.0 ^ M >-22.0 


1.0 


3834 


1.159±0.009 


0.232±0.001 


-6.738±0.195 


0.274 


-22.0 ^ M >-23.0 


1.0 


1824 


1.220±0.012 


0.244±0.001 


-7.062±0.081 


0.259 


-23.0 ^ M >-24.0 


1.0 


144 


1.330±0.055 


0.266±0.001 


-7.719±0.371 


0.228 








BCESiog(r^) 








-18.5 > M >-19.0 


0.5 


26 


1.267±0.111 


0.253±0.001 


-7.630±0.691 


0.062 


-19.0 ^ M >-20.0 


1.0 


492 


1.591±0.079 


0.318±0.001 


-9.562±0.492 


0.083 


-20.0 ^ M >-21.0 


1.0 


2341 


1.504±0.031 


0.301±0.001 


-9.000±0.199 


0.083 


-21.0 ^ M >-22.0 


1.0 


3834 


1.458±0.024 


0.292±0.001 


-8.673±0.158 


0.075 


-22.0 ^ M >-23.0 


1.0 


1824 


1.448±0.032 


0.290±0.001 


-8.569±0.212 


0.069 


-23.0 ^ M >-24.0 


1.0 


144 


1.480±0.122 


0.296±0.001 


-8.726±0.826 


0.062 








BCES<,i>^ 








-18.5 ^ M >-19.0 


0.5 


26 


1.492±0.019 


0.244±0.016 


-6.900±0.326 


0.312 


-19.0 ^ M >-20.0 


1.0 


492 


1.771±0.090 


0.266±0.001 


-7.245±0.141 


0.435 


-20.0 ^ M >-21.0 


1.0 


2341 


1.731±0.044 


0.263±0.001 


-7.107±0.069 


0.429 


-21.0 ^ M >-22.0 


1.0 


3834 


1.630±0.036 


0.255±0.001 


-6.864±0.059 


0.395 


-22.0 ^ M >-23.0 


1.0 


1824 


1.586±0.048 


0.252±0.001 


-6.707±0.081 


0.363 


-23.0 ^ M >-24.0 


1.0 


144 


1.613±0.163 


0.254±0.013 


-6.659±0.276 


0.322 
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respect to the log(re), M and Xfp variables shows a corre- 
lation. In particular for the magnitude case the correlation is 
very clear. Therefore, the residue distribution confirms the 
fact that different luminosity galaxies are found in consec- 
utive planes characterised by different coefficients. But why 
do the values of the coefficients of the FP for different lumi- 
nosity galaxies change? Section 3.5 will deal in greater detail 
with this. 

3.4 Previous works dealing with the behaviour of 
the FP with respect to magnitude 

In the literature there are several papers which study the 
dependence of the FP with magnitude, for example Jor- 
gensen et al. (1996) and Bernardi et al. (2003c); utilising 
the residues of the FP they find that there is no dependence 
of the FP with magnitude. However, the typo of analysis 
undertaken by them might be masking the dependence on 
the magnitude range which we are reporting in this paper. 
Given that in the Bernardi et al. (2003c) paper they use the 
exact same sample of ETGs from the SDSS as we do here, 
in what follows we shall analyse in more detail their results. 

Among the most important conclusions of Bernardi et 
al. (2003c) related to our work, we find that the orthogo- 
nal residues of the FP (when they take as reference the or- 
thogonal fit coefficients) with respect to the log (re), < fj. >e 
,log(a), M and Xfp variables do not show an appreciable 
correlation, except for the log(cr) and Xpp cases. On the 
other hand, when they analyse the residues taking the di- 
rect fit FP coefficients and using as independent variable 
log(re), they mention that no correlations are found. They 
interpret the correlation found for the orthogonal fit as a pro- 
jection effect which affects this type (orthogonal) of fits but 
not the direct fits (if log(re) is taken as dependent variable). 
However in figure 5 we show the results of the orthogonal 
residues (g* filter) taking the data from the direct fit (log (re) 
is the independent variable) from the table 3 of Bernardi 
et al. (2003c) and it can be clearly noticed that there is a 
correlation of the residues with respect to the log(re) and 
magnitude. Particularly for the magnitude case the correla- 
tion is perfectly clear and apparent. On the other hand, the 
data from table 4 show that the change in the coefficients 
of the FP occurs no matter which variable is utilised as in- 
dependent variable in the fitting process. So we are able to 
ensure that the fitting method and the independent variable 
utilised for the fits are not responsible for the dependence 
of the FP on the magnitude range which we are reporting 
in this paper. 

The distribution of the residues for the different vari- 
ables presented in figure 3 has been obtained in a similar 
manner as that followed by Bernardi et al. (2003c), that is 
to say, taking as reference the FP coefficients for the total 
sample and taking into consideration the different sources 
of bias. When we compare these distributions with those 
presented in Figs. 2 and 3 of Bernardi et al. (2003c), we no- 
tice that the residues distribution is similar, that is, there 
are no appreciable correlations between the values of the 
residues and the different variables which we analyse (except 
for the log(o') and Xfp cases). However, as wo mention in 
Section 3.3, when we take as reference a representative plane 
of the sample in one-magnitude wide intervals, for example 
the plane at the middle part of the distribution of galax- 



ies (-21.0 ^ M >-22.0) (see figure 4), the residues clearly 
show a correlation with respect to the log(re), M and Xfp 
variables, in particular, the residues show a clear correlation 
with respect to magnitude. From all of this we can say with 
certainty that in the Bernardi et al. (2003c) analysis the de- 
pendence of the FP with the range of magnitude has been 
masked due to the way in which these authors performed 
the analysis. 



3.5 Comparison of the behaviour of the FP and 
the KR with respect to the magnitude range 

Nigoche-Netro (2007) and Nigoche-Netro et al. (2007, 2008) 
demonstrated, both with observational data as well as with 
numerical simulations that the coefficients and the intrinsic 
dispersion of the KR depend on the width and brightness of 
the magnitude range. This dependence is caused by a "geo- 
metrical effect" due to the fact that the distribution of the 
ETGs on the log(re) - < >e plane depends on luminosity 
(different brightness ETGs are distributed in parallel lines 
on the log(re) - < jj, >e plane) and that the geometric shape 
of the distribution of the ETGs on this plane is not sym- 
metrical. Given this fact, and since the KR is a projection 
of the FP, the dependence of the values of the FP coeffi- 
cients on the width and brightness of the magnitude range 
might be caused by a geometrical effect due to the form of 
the distribution of the ETGs on the space defined by the 
variables log(re), < H >e,log(cr) and not by intrinsic physi- 
cal properties of the ETGs. In fact this is the case given that 
in Section 3.3 we have seen that the distribution of ETGs 
in the log(re),< fi >e,log((T) space depends on luminosity, 
besides in Section 3.2 and 3.3 we have established that the 
values of the coefficients and of the dispersion of the FP 
change systematically as we consider brighter galaxies. This 
implies that the geometric shape of the galaxy distribution 
changes as brighter galaxies are included. This may easily 
be checked from figure 3 where one can see that as brighter 
galaxies are considered the shape of the galaxy distribution 
changes systematically. 

From the aforementioned, we may infer that the FP 
coefficients and its intrinsic dispersion depend on the width 
and brightness of the magnitude range. This dependence 
is caused by a geometrical effect due to the fact that the 
distribution of the ETGs on the log(re) , < /i >e, log((7) space 
depends on luminosity, in other words, the geometric form 
of the distribution of ETGs on the log(re),< n >e,log(cr) 
space changes systematically as we consider brighter ETGs, 
and so the values of the FP coefficients change too, because 
the fitting of a plane to a set of data docs not give the same 
result for data distributed with a cubic form as for data 
distributed with a piramidal form or, for that matter, with 
another geometrical form. In this sense, any other systematic 
restrictions imposed on a sample of ETGs, such as brightness 
cuts, effective radius cuts or velocity dispersion cuts will 
cause changes in the geometric form of the distribution of 
ETGs. The more pronounced these changes are made, the 
more pronounced will be the changes in the values of the 
FP coefficients. This could be corroborated with the data 
in figure 1, where we show that the coefficients of the FP 
for the homogeneous sample are larger than those for the 
complete sample. This behaviour is due to the cut in log(cr) 
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Figure 3. Orthogonal residues for the FP vs. log(re),M, < fi >e,log(o'), and Xpp. The FP coefficients we use correspond to the 
BCESsis fit for all the galaxies in the SDSS sample (g* filter). Each symbol and colour represent a one-magnitude wide interval. 




Figure 4. Orthogonal residues for the FP vs. log{re), M, < fi >e,log((T), and Xpp. The FP coefficients we use correspond to the 
BCESpis fit for the middle part (-21.0 M >-22.0) of the SDSS galaxy distribution (g* filter). Each symbol and colour represent a 
one-magnitude wide interval. 
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Figure 5. Orthogonal residues for tlie FP vs. log(re),M, < ^ >e,log(fr), and Xpp. The FP coefficients correspond to the direct fit 
(log(re) independent variable) for Table 3 from Bernardi et al. (2003c) (complete sample of the SDSS in filter g*). Each symbol and 
colour represent a one-magnitude wide interval. 



and not to the fact that the FP for this galaxy subsample 
might be intrinsically different. 

In view of what we present above, the traditional FP 
relation (equation 1) could be written as follows: 

log (re) = a'log(cro) + b' < >e + c (3) 

where a' , b' and c' depend on the width and brightness 
of the magnitude range. 

It is very important to be careful with equations 1 and 3 
because, in order to be consistent, the first one should cover 
the whole range of magnitudes and the second should corre- 
spond to a particular magnitude range. In other words, the 
traditional FP relation will only produce the real coefficients 
when complete samples, as far as magnitude is concerned, 
are considered. 

Here it should be mentioned that in the literature there 
are extensive studies of the FP for different galaxy samples. 
The values of the coefficients of the FP which are reported 
therein show important differences among the various au- 
thors. For example D'Onofrio et al. (2006) have made a 
study of several papers that have dealt with the FP. From 
the D'Onofrio et al. (2006) data we find that the maximum 
difference of the values of the a coefficient estimated by dif- 
ferent authors is of the order 31%, for the b coefficient the 
maximum difference amounts to 16% and for the c coeffi- 
cient there is no information given. In our paper, when we 
compare the values of the coefficients for galaxy samples in 
narrow magnitude intervals with the values obtained from 
galaxy samples in wide magnitude intervals, we find that the 
maximum difference in the a coefficient is approximately 
of the order 33%, for the b coefficient it is approximately 
10% and for the c coefficient it is 20%. Given the similarity 



between the differences in the coefficients in our work and 
those obtained using the data from D'Onofrio et al. (2006), 
we assert that an important reason for the differences found 
among the various authors may be caused by the different 
magnitude ranges within which the galaxies are contained, 
on top of the differences produced by the fitting method 
utilised and the independent variable used in this fitting 
method. On the other hand given that the differences found 
for the values of the coefficients a and b in our work are 
similar to those reported in the literature, we can safely say 
that the differences for the values of the c coefficient from 
the literature will also be similar to the differences in the c 
values reported here, approximately 20%. 

An important application of the FP is that it may be 
used to perform large scale distance estimates, however, the 
differences between the distances calculated utilising the val- 
ues of the FP coefficients from different authors reported in 
D'Onofrio et al. (2006), may be significant. For example, if 
we assume that the FP is universal, that is to say, that the 
values of the coefficients a and b are consistent with a unique 
value, then the distances would be given as function of the 
differences in the zero point (see Blakeslee et al. 2002), then, 
if there is a 10% error in the zero point this would produce 
a 50% error in the distance determination. This is the error 
which could be made if the data reported in D'Onofrio et 
al. (2006) were used and it could also be the error made if 
the dependence on the magnitude range were not consid- 
ered. Here we must stress that these distance uncertainties 
have been obtained without considering in the calculation 
the errors in coefficients a and b, so the total error could be 
larger than 50% both if we use the FP data from the differ- 
ent authors reported by D'Onofrio et al. (2006) as well as if 
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we do not take into consideration the dependence of the FP 
coefficients on the magnitude range. 

Finally, due to the change of the traditional FP coeffi- 
cients with respect to the width and brightness of the mag- 
nitude range, it is important to consider that when com- 
parisons between different ETGs samples are performed, it 
must be clearly established what is the magnitude range 
within which the ETGs are distributed. If this is not done 
the differences which might be found may be misinterpreted. 

4 CONCLUSIONS 

In this paper we make a compilation of photometric and 

spectroscopic parameters for 3 samples of ETGs from the lit- 
erature, which include a total of approximately 8800 galax- 
ies and cover a relatively ample magnitude range (< AM > 
^ 5.5 mag). With this information we perform an analysis 
of the behaviour of the coefficients of the FP with respect 
to the width and brightness of the magnitude range. The 
results we obtain are as follows: 

• We find that when a growing number of brighter galax- 
ies is included in the calculations of the intrinsic dispersion 
and of the coefficients of FP (increasing magnitude intervals) 
or if we consider galaxy samples in progressively brighter 
fixed-width magnitude intervals (narrow magnitude inter- 
vals), their values change and these changes are larger than 
the associated errors for most of the cases. 

• This fact made us think that there might be a depen- 
dence both of the values of the coefficients of the FP rela- 
tion and of the intrinsic dispersion of this relation with the 
width and brightness of the magnitude range. The existence 
of this dependence is confirmed with the application of a 
non-parametric test to the data (run-test). 

• Since the coefficients of the FP relation depend on the 
width and brightness of the magnitude range, we establish 
that this relation could be written as follows: 

log (re) = o' log ((To) + h' <M>e + c' (4) 

According to this, the traditional FP relation (equation 
1) will produce the real coefficients only when we consider 
complete samples, as far as magnitude is concerned. 

• We find that the behaviour of the FP relation described 
above is similar to the bahaviour of the KR reported in 
Nigoche-Netro (2007) and Nigoche-Netro et al. (2007, 2008), 
that is, the dependence of the FP on the width and bright- 
ness of the magnitude range is caused by a geometrical effect 
due to the fact that the distribution of the galaxies in the 
log(re), < /U >e,log((T) space depends on luminosity. 

• We find that any systematic restrictions imposed on a 
sample of ETGs, such as brightness cuts, effective radius 
cuts or velocity dispersion cuts will cause changes in the 
geometric form of the distribution of ETGs on the log(re), < 
fj, >e,log(cr) space. The more pronounced these changes are 
made, the more pronounced will be the changes in the values 
of the FP coefficients. 

• Finally, due to the change of the FP coefficients with 
respect to the width and brightness of the magnitude range, 
it is important to remember that when comparisons between 
galaxy samples are performed, these must be done only for 
equivalent magnitude intervals, otherwise the results may be 
misinterpreted. 



ACKNOWLEDGMENTS 

We would like to thank Consejo Nacional de Ciencia y 
Tecnologia (Mexico) for a PhD fellowship number 132526, 
Ministerio Espaiiol de Educacion y Ciencia for grant 
PNAYA2006, Institute de Matematicas y Ffsica Fundamen- 
tal (CSIC, Espana) and Instituto de Astronomia (UNAM, 
Mexico) for all the facilities provided for the realisation of 
this project. We would also like to record our obligation to 
Prof. Mariano Moles for help with the presentation of this 
paper. Last but not least, we would like to acknowledge the 
many suggestions provided by an anonymous referee which 
greatly improved the presentation of this paper. 



REFERENCES 

Akritas, M. G., & Bershady, M. A. 1996, ApJ, 470, 706 
Bondat J. S., & Picrsol A. G. 1966. Measurement and Anal- 
ysis of Random Data. Ed. Wiley J. & Sons. New York. 
Bender, R., Burstein, D., & Faber S. M. 1992, ApJ, 399. 
462 

Bernardi, M., et al. 2003a, AJ, 125, 1817 

Bernardi, M., et al. 2003b, AJ, 125, 1849 

Bernardi, M., et al. 2003c, AJ, 125, 1881 

Bernardi, M., et al. 2007a, A J, 133, 1741 

Bernardi, M., et al. 2007b, A J, 133, 1954 

Blakeslee, J. P., Lucey, J. R., Tonry, J. L., Hudson, M. J., 

Narayanan, V. K., & Barris, B. J. 2002, MNRAS, 330, 443 

Colberg, J. M., et al. 2000, MNRAS, 319, 209 

Djorgovsky, S., & Davies, M. 1987, ApJ, 313, 59 

D'Onofrio, M., Valentinuzzi, T., Secco, L., Caimmi, R., & 

Bindoni, D. 2006, NewAR, 50, 447 

Dressier, A., Lynden-Bell, D., Burstein, D., Davies, R. L., 
Faber, S. M., Terlevich, R., & Wegner, G. 1987, ApJ, 313, 
42 

Efron, B. 1979, Ann. Statist., 7, 1 

Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shi- 
masaku, K., & Schneider D. P. 1996, AJ, 111, 4 
Hamabe, M., & Kormendy, J. 1987, In: Structure and Dy- 
namics of Elliptical Galaxies, lAU Symp., No 127, p. 379, 
ed. de Zeeuw, T., Reidel, Dordrecht 

Hoessel, J. G., Oegerle, W. R., & Schneider, D. P. 1987, 
AJ, 94, 1111 

Isobe, T., Feigelson, E. D., Akritas, M. G., & Babu, G. J. 
1990, ApJ, 364, 104 

Jorgensen, I., Franx, M., & Kaergaard, p. 1996, MNRAS, 
280, 167 

Jorgensen, I., Franx, M., Hjorth, J., & van Dokkum, P. G. 

1999, MNRAS, 308, 833 

Kelson, D.D., van Dokkum, P. G., Franx, M., lUingworth, 
G. D., & Fabricant, D. 1997, ApJ, 478, L13 
Kjasrgaard, P., Jorgensen, I., & and Moles, M. 1993, ApJ, 
418, 617 

Kormendy, J., 1977, ApJ, 218, 333 

La Barbera, F., Busarello, G., & Capaccioh, M. 2000, A&A, 
362, 851 

La Barbera, F., Busarello, G., Merluzzi, P., Massarotti, M., 

& Capaccioh M. 2003, ApJ, 595, 127 

Milvang- Jensen, B. 1997, Master's Thesis. University of 

Copenhagen 

Nigoche-Netro, A., Moles, M., Ruelas-Mayorga, A., Franco- 
Balderas, A., & Kjaergaard, P. 2007, A&A, 472, 773. 



The Fundamental Plane for ETGs. Dependence on the magnitude range 11 



Nigoche-Netro, A. 2007, PhD Thesis. Universidad Com- 
plutense de Madrid 

Nigoche-Netro, A., Ruelas-Mayorga, A., & Franco- 

Balderas, A 2008, A&A, Accepted (arXiv:0805.0961 ) 

Prugniel, P., & Siemen, F. 1997, A&A, 321, 111 

Reda, F. M., Forbes, D. A., Beasley, M., O'Sullivan, E. J., 

& Goudfrooij, P. 2004, MNRAS, 354, 851 

Sandage, A., & Peremulter, J-M. 1991, ApJ, 370, 455 

Sandage, A., & Lubin, L. M. 2001, ApJ, 121, 2271 

Varela, J. 2004, PhD Thesis. Universidad Complutense de 

Madrid 



APPENDIX A: ANALYSIS OF THE 
ROBUSTNESS OF THE RUN TEST USING THE 
BOOTSTRAP METHODOLOGY 

Al Run test robustness analysis for a set of 
general data 

We have run a number of experiments with the run test. 
This is the test which we use to establish the presence of an 
underlying trend in the behaviour of the coefficients of the 
FP. The experiments consist in: 

• Generating a list of 200 random numbers between and 
1. A test of their distribution with the run test reveals that 
there is absolutely no underlying trend. 

• Following the Bootstrap method, we obtain 1000 resam- 
ples from the original sample and test their distribution with 
the run test. The character of their distributions changes 
from not having an underlying trend to having one in less 
than 20 cases in 1000 at a confidence level of 99%. 

• Artificially introducing a monotonic trend. 

• Testing the new distributions with the run test reveals, 
in all the cases, the presence of an underlying trend at a 
confidence level of 99%. 

• Now we introduce randomly an error which may be pos- 
itive or negative and we test the distributions again. The 
errors introduced were ±0.1, ±0.2, ±0.3, ±0.4 and ±0.5 or 
errors of the order 10%, 20%, 30%, 40% and 50%. 

• For all the cases except for that with 50% error, the 
run test continued to indicate the presence of an underlying 
trend at a confidence level of 99%. 

• Again following the Bootstrap method, we obtain 1000 
resamples from the original sample with random errors. 

• The run test reveals that in these thousand resamples 
with errors, the character of the distributions changes from 
having an underlying trend to not having one in less than 
20 cases out of 1000 at a confidence level of 99%. 

• The run test reveals that if the number of data for the 
original sample decreases and the error in the data increases 
then the number of resamples that change their distribution 
from having an underlying trend to not having one grows 
and the level of confidence at which we may reject the null 
hypothesis becomes smaller. 

We conclude from this analysis that the run test predicts 
reliably the presence of an underlying trend in a set of data, 
even if these data have errors. Only when the errors are very 
large (>50%), and the amount of data is small the run test is 
unable, in some cases, to detect the underlying trend present 
in the data. 



Table Al. Run test for the evaluation of the 
FP coefficients (increasing magnitude intervals) 
from the different samples of galaxies. One thou- 
sand samples were extracted from the original 
data following the bootstrap methodology. In 
addition the data were combined in a random 
way with Gaussian errors for which the max- 
imum errors corresponded to the values cited 
in Table 2. The null hypothesis establishes that 
there is no underlying trend in the data. The 
percentages refer to the fraction of extracted 
samples for which an underlying trend is de- 
tected at a confidence level of 90%. 



Sample 


a 


b 


c 


Total SDSS (g* filter) 


94% 


94% 


96% 


Total SDSS (r* filter) 


94% 


93% 


96% 


Total SDSS (i* filter) 


95% 


93% 


95% 


Total SDSS (z* filter) 


95% 


95% 


95% 


Coma (Gunn r filter) 


82% 


82% 


80% 


Hydra ( Gunn r filter) 


82% 


95% 


80% 



A2 Run test robustness analysis for the FP data 

We ran a series of experiments on the values of the coef- 
ficients of the FP using the bootstrap methodology. One 
thousand samples were extracted from the original data (see 
Table 2). We also added to the data Gaussian errors. The 
maximum allowed error for each datum corresponds to the 
error for the real data. In Table Al we show the results of 
the run test where the percentage given refers to the fraction 
of extracted samples in which an underlying trend has been 
detected at a level of confidence of 90%. 

As can be seen in Table Al, when the maximum allowed 
error in the data of the extracted samples is of the same size 
as those errors affecting the real samples, we find that on the 
average 91% of the extracted samples show an underlying 
trend at a confidence level of 90%. We may again affirm 
that the run test is efficient in finding underlying trends 
even when the errors in the data have a significant value. 

This paper has been typeset from a TgX/ file prepared 

by the author. 



